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We obtain a class of adiabatic solutions of Dirac equation for the charged massless relativistic 
quasi-particles that arise from the low-energy excitations p_i in a 2D graphene sheet, interacting with 
an electromagnetic field. The analytic solutions obtained are useful for non-perturbative investigation 
of processes in intense laser fields. As a first example we employ them to predict copious emissions 
of high order harmonics by THz lasers interacting with the occupied states of graphinos in the 
vicinity of the degenerate Dirac points. The relative intensity of the emitted harmonics is seen to 
decrease by only about two orders of magnitude from the 3rd to the 81st harmonic order, and is 
characterized by two phenomena of "revival" and "plateau" formation in the middle and the far 
end of the spectrum calculated. A preliminary comparison is made for harmonic emission from 2D 
graphene that reveals a qualitatively different spectrum in the latter case showing a sharp cutoff at 
an order n cu toff = ^j- where lob = eF^b/h is the so-called Bloch frequency. 

PACS numbers: 78.67.Wj, 78.70.-g, 78.47.jh,03.65.Xp 



INTRODUCTION 



Graphene is a two-dimensional single layer of carbon 
atoms in a hexagonal honeycomb arrangement (Fig. 1 
(a)). Since its recent controlled synthesis in the labora- 
tory graphene has emerged as an object of much current 
research interest [3] with potentials for wide range of 
applications in various quantum technologies. A remark- 
able property of a 2D graphene membrane is that its 
low energy excitations - near the so-called Dirac points 
(K, K') at the corners of the graphene Brillouin zone (cf. 
Fig. 1 (b)) - are relativistic fermions that are charged 
like the electrons but massless like the neutrinos. They 
move with a constant velocity vf that makes them fully 
relativistic at a much lower velocity than the velocity of 
light c (vp s» c/300). Also, their anti-particle counter- 
parts are oppositely charged massless holes with opposite 
helicity [1]. 

The stationary eigenvalues and the spinor eigenfunc- 
tions of the Dirac equation of free graphinos as well as the 
wavefunctions in the presence of static electric and mag- 
netic fields, and in crossed electric plus magnetic fields 
have been obtained recently (see e.g. review[3j and [HH]). 
They have been used to investigte some of the most re- 
markable properties of graphinos, such as chiral scatter- 
ing, confinement, zitterbewegung, and anomalous integer 
quantum Hall effect - among others - that have been ob- 
served [2j |3] . In this paper we consider the adiabatic so- 
lutions of the time- dependent Dirac equation of graphinos 
interacting with an electromagnetic field of a given pulse 
shape, peak field strength, central frequency, and polar- 
ization. As a first application of the solution obtained, we 
derive the Dirac expectation value of the graphino cur- 
rent in the presence of the field and show that copious 
emission of high harmonic radiations can occur from the 
initially occupied states of graphino pseudo-particles near 




Figure 1: Structure of graphene: (a) 2D hexagonal honey- 
comb lattice, ai = ax + by, a,2 = ax — by, a — |ao, b — -^ao; 
(b) Brillouin zone showing the Dirac points K(-^, 3 ^ a ) 
and K' 
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the degenarate Dirac points, by an intense THz laser. 

DIRAC EQUATION OF A GRAPHINO 
INTERACTING WITH A TRANSVERSE 
ELECTROMAGNETIC FIELD 

The ID Dirac Hamiltonian of a graphino is given by 
H = v f <t ■ p op (1) 
where a = (cr x , o y ) are 2x2 Pauli matrices 



1\ _ (0 -i 

1 j ' Uy ~ [i 



(2) 
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v F is the graphino velocity, and p op = —ih-^, r = (x, y). 
We consider a general plane wave electromagnetic field, 
assumed to be elliptically polarized in the graphene 
plane, to be represented by the vector potential 

A W = Aoirj) (e^ cos (£/2) cos 7] Cy sin (£/2) sin fj) , 
rj = (wt — kz) (3) 

where Ao(rf) is the pulse envelope, £ = [0, ±7r/2] is the el- 
lipticity parameter with left (+) or right (— ) handedness, 
and e x ,e y are unit polarization vectors. The minimal- 
coupling prescription yields the Dirac equation of the 
"graphino + electromagnetic field" : 



The unperturbed graphino 



(4) 



It will be useful first to briefly consider the eigen- 
solutions [3] of the field-free graphino Dirac equation gov- 
erned by the Hamiltonian ([!]), 



Eip p (r) = v F a ■ p^p(r) 



(5) 



where E is the energy and p — hq is the quasi-momcntum 
measured from the Dirac point K (or K'). One substi- 
tutes Vp( r ) — es p r u p , where 



(6) 



is a two component spinor, in Eq. ^ to get 
Eu p = vpo ■ pu p 



v F p 











(7) 



where (f> p is the polar angle of the 2D momentum p. This 
gives a 2 x 2 matrix equation for the components a and 
b of the spinor u p 

e-'M^, (8) 



E 



b) =VFP \e^ 







This can be readily solved to obtain the two eigenvalues 

E {±) = ±v F p 
and the corresponding normalized eigenspinors, 

M ( ±) _ JL 

P ~ V2 \±e 1 ^/ 2 J ■ 
Thus the normalized eigenstates of Eq. ^ are 

I 



(9) 
(10) 



V£)(r) = 



e iP' r u ( s ) 
^ p 



(11) 



with eigenvalues E^ s '(p) — sv F p, where A is the quanti- 
zation area, and s = ±1 are the two components of the 
graphino psedo-spin, corresponding to the limiting case 
of the upper and the lower band of graphene, respectively. 



Adiabatic solutions of graphino Dirac equation in a 
laser held 



We may now proceed to obtain a class of useful solu- 
tions of the time-dependent Dirac equation Q for the in- 
teracting system "graphino + laser field" . We first make 
the following adiabatic ansatz for the wavefunction (e.g. 

i) 

¥ P (r, t) = -^e-T>f E ^ dt ' e^ r u p (t) (12) 
v A 



where E(t) and u p (t) are time dependent unknowns to 
be determined. To this end we substitute Q in Eq. Q 
and obtain 



d 

E(t)u p {t) + ih—u p (t) 



= (v F a-(p-^A(r,)))u p (t) (13) 

Next, projection from the left with the hermitian adjoint 
spinor u p (t) gives 

E(t)ul(t)u p (t) + ihul(t)u p (t) 
= «t( t ) (v F o ■ (p - - c Hv)i) «p(*) (14) 

where we have denoted u p (t) = J^w p (t). Comparing the 
last term in ( 14 ) with the unperturbed equations Q and 
( 10 1 we make the ansatz 



«p(*) = « P (t) = ^ \ se i<t>H^) ' < 15 ) 

where we have defined p(t) = (p — |A(rj)),p(t) = |p(i)|, 
and tan</) p ( t ) = p y (t)/p x (t). It is easily seen by direct 
calculation that, 



(«)t (*') _ V 



(16) 



By calculations similar to that shown above for the un- 
perturbed case but with the time-dependent spinors de- 



fined by Eq. ( 15 ) we also find 



(v F a ■ (p - e c A(rj))) u<*>(*) = sv F \p(t)\u^(t)(l7) 
Moreover, since 



«W(t) 



L ^" i0p(t)/2 (-|^(t)) 



72 V se^w/ 2 (|0 p(t) ) 



(18) 



and 



uWt(t) = ^( e ^ P (*)/ 2 , se^W 2 ), (19) 



therefore the inner product 



1 



uW(t)u p s >(t) = -^-i^+i^t)) = (20) 
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Thus, substituting Eqs. (16 1, (17 1 and (20) in Eq. (14) and from its hermitian conjugate 



it is seen that Eq. ([14]) is fully satisfied provided 

E(t) = sv F \p(t)\ (21) 

Hence, finally, we obtain the desired adiabatic solution 
of the time-dependent Dirac equation Q 

<?W(r,i) = JL e -farr\ptf)\tf e iP-' u W (22) 
v A 

with s = ±1. They fulfil the biorthonormal conditions 

J ^ 3t (r,t)*M(r,f)d 2 r = 5(p-p')6 a , s , 
X;^*W(r',t)*Wt(r,i) = J(r-r') (23) 

p s— ± 

and thus form a complete set of functions that can be 
used to represent an adiabatic state |$(t)) satisfying any 
given initial condition |<j>(io)) at a time to, in the form 

= E E<*p } (*o)i*(*o)) |* P s) w) ( 24 ) 

p S— ± 



HIGH HARMONIC EMISSION FROM 
GRAPHINO 

As a first application of the solution obtained above, 
we consider generation of higher order harmonics from 
graphinos near the Dirac points by intense laser fields. 
Emission of high harmonic of radiation might be ex- 
pected in analogy with the earlier theoretical investiga- 
tions within the elaborate Floquet-Bloch theory for crys- 
tal surfaces and graphene [IrlS]) or from calculations with 
two-state models of semiconductors [TU] , as well as from 
experimental observations of higher order harmonic emis- 
sions from surfaces and bulks [TTHT3"] by high intensity 
infrared laser fields. In view of the zero-mass of graphi- 
nos we may expect their strong effective coupling with 
low frequency laser fields and consequent copious emis- 
sion of high harmonics, even by moderately intense lasers. 
To show this theoretically we next obtain the quantum 
expectation value of associated current operator in an 
electromagnetic field and calculate the relative harmonic 
intensities from a state occupied in the vicinity of the 
Dirac points, for the case of a THz laser. 



Graphino current driven by an electromagnetic field 

From Eq.([5]) we have 
d 



dt 



*p(t) - -UF<7- W p (t) 



|U P (i) = -v F a-V%(t) 
= -v F \7 ■ ^l(t)a 



(26) 



since ^ = a. Multiplying (25) with # P (i) and (26 1 with 
^p(t) and adding them we obtain 



0_ 
dt 



(t)* P (t)) = - VF V ■ (*;(i) ff * p (f)) (27) 



or, the graphino continuity equation: 
d 



dt 



p(f)+V-j(t) = 



(28) 



where 



p(t) = * P («)*p(i) 



(*)t W 
u P (t) u P (t) 



is the probability density, and 



j(t) = v F yUt)aV p (t) 



0)t 0) 

p(t) au p(t) 



= v F u 



(29) 



(30) 



is the probability current density. We assume that the 
lower band (s= -1) is fully occupied and therefore in an 
adiabatic laser field harmonic emission is due essentially 
to the coherently driven graphino current in the upper 
band (s=+l), arising from the initially occupied states 
near the degenerate Dirac points. For a laser incident 
perpendicular to the graphene plane and linearly polar- 



ized along the x-axis, we substitute the spinors ( 15 ) and 



(19) in Eq. (30) and readily evaluate the quantum ex- 



pectation value of the graphino current operator to get: 



4 s) (P,i) 



-ej x (t) 

-esv F cos 4> p {t) 

Px - -A x (ri) 



-esvp 



^( Px -IA.M) 2 +Pl 



(31) 



Similarly, for the laser polarization along the y-axis we 
obtain 

Jy s) (p,t) = -esv F sin <£ p (t) 

P y - - c A y (r)) 
= ~ esv F—p^ , pa i \ \ 9 (32) 

VpI + (Py~ -cMv)) 2 

and for the more general case of an elliptically polarized 
laser field (cf. EqM) we get [H]: 



(25) 



J (s) (p,£) = - esv F (e x cos c/> p (t) + e y sin </> p (t)) 

e x (Px - - c A x {rj)) + e y (p y - %A y (rj)) 

^-^(i)) 2 + (p»-fW 

(33) 



= —esvp- 
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Figure 2: High order harmonic emission (shown up to the 
13th harmonic) by THz laser interaction with "graphinos" 
near the Dirac point K(-^r-,^w — ); illustrative calculation 



for 



(K x 



•0 ' 3\/3ao ; ' 

0.01, q y b 



(Ky 



ky)b 



0.01, 



laser polarization along y-axis, frequency 3 THz, focused laser 
intensity I = 50 MW/cm 2 . Note the presence of the "revivals" 
and "plateaus" in the middle and the end of distribution. 



where A x (rj) = A (rj) cos (£/2) cos?7, and A y {rf) — 
Aq(t]) sin (£/2) sin rj. The intensity of the harmonic emis- 
sion at the nth harmonic frequency £1 = nuj depends es- 
sentially on the absolute square of the Fourier transform 
of the graphino current: 



S (s \nio) = 



(nuj) 



f e -«»Jj(s)(p ;?7 ) rf?? 

J —TV 



(34) 



where n = 1, 2, 3, • • • . In Fig. 2 we present the re- 
sults of calculations for the relative intensity of harmonics 
emitted from an occupied state in the neighborhood of 
the Dirac point K, by a 3 THz laser with a focused in- 
tensity of 50 MW/cm 2 , that is polarized linearly along 
the y-axis (cf. Fig. 1) and is directed perpendicular 
to the graphene plane. It should be noted that exactly 



at a Dirac point (e.g. K = {-, 



-)) the driven 



k 3a ' 3V3a ' 

graphino current in the present case becomes a constant 
(cf . Eq. ( 32 1 ) and therefore no harmonics are pro- 
duced by graphinos from this point. Nevertheless, as 
expected above high order harmonic emission can occur 
copiously from the initially occupied states in the upper 
band (s=+l) in the vicinity of the Dirac points. Fig. 2 
shows an example of harmonic emission from a neighbor- 
ing point (q x a = 0.01, q y b = 0.01), exhibiting highly non- 
perturbative relative intensity [TS] that is seen to decrease 
only about two orders of magnitude over some 81 orders 
of the harmonics. Moreover, it can seen clearly from the 
figure that the distribution of relative intensity is charac- 
terized by "revivals" and multiple "plateaus" within the 
range shown in Fig. 2 [15] , 



Harmonic generation from graphene band 

The full dispersion relation of graphene is very differ- 
ent from that of the graphino quasi-particles. It should 
be therefore interesting to see the difference between the 
harmonic generation spectra in the two cases. We discuss 
below the results of our preliminary investigation of the 
difference in a particular case. To this end we consider 
the energy band of a graphene sheet in the tight-binding 
approximation (see, e.g. [3]). In this approximation the 
upper and the lower bands are given 

E {s \p) = st{3 + 2 cos (2k y b) + 4cos(k x a) cos(k v b)} 5 
= st{l + 4cos (k x a) cos (k y b) + (2cos (k y b)) 2 }^ 

(35) 

The semiclassical band current can be obtained by dif- 
ferentiating the band energy with respect to the band 
momentum and applying the minimal coupling prescrip- 
tion. Thus for the band current induced by an elliptically 
polarized electromagnetic field propagating perpendicu- 
lar to the graphene plane we find, 



J [s) {p,t) 



Jy S \p,t) 



st(—2a) sin (k x (r))a) cos (k y (r))b) 
x{l + 4cos (k x (r])a) cos (k y (rj)b) 
+ (2cos(fc,,(77)&)) 2 }-5 

(-26){cos (k x (r])a) sin {k y {rf)b) 
+ cos2(k y (r 1 )b)} 

x{l + 4cos (k x (r])a) cos (k y (r])b) 
+ (2cos(fc,(77)&)) 2 )}-5 



(36) 



(37) 



where, 



k (v) = (k x (r)),k y (r))) 



k x (r)) 



he 



A x (0 cos (r)) 



k y (r)) = k y - —Ay((,) sin (ri) 



(38) 
(39) 



Let the laser be polarized along the y-axis and assuming 
that the initial states along k y and k x = contribute the 
most to the driven current in this case we get: 

3 y s) {Py,t) = (-2stb) sin (k y (r,)b) 

= (-2stb) sin (k y - p sin 77) (40) 

where, p = and Fq = ^§ SI - Using the Jacoby- Anger 
relation 



e ipmn V= J2 J n {p)e mV 



(41) 
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Figure 3: High order harmonic emission (shown up to the 
46th harmonic) by aTHz laser interaction with a 2D graphene 
layer near the Dirac point K(^-, —7= — ); illustrative calcu- 
lation for q x a = (K x — k x )a — 0, q y b = (K y — k y )b = 0.01, 
laser polarization along y-axis, frequency 3 THz, focused laser 
intensity I = 500 GW/cm 2 . Note the presence of a cut- 
off order n cu toff ~ 20 after which the relative harmonic 
intensity drops exponentially, consistent with the formula 
ncutoff « J = 19.9. 



where e ^j^- is the so-called Bloch frequency. Fi- 
nally, in Fig. 3 we show the relative intensity distribution 
of the harmonics for the case of a 3 THz laser, polarized 
along the y-axis, at an intensity I = 500 GW/cm 2 , near 
the Dirac point (q x a = 7 q y b = 0.01). As can be seen, 
now the harmonic spectrum extends significantly only 
up to about n cuto ff ~ 20 that confirms the prediction 
of the cutoff- formula n cuto ff ~ — 19.9; beyond this 
value it indeed falls of exponentially. This behavior is 
qualitatively different from that of the free graphinos dis- 
cussed above (cf. Fig. 2), as well as from the case of free 
atoms and molecules for which a cutoff order exists but is 
given differently by the formula: n cuto ff ~ I p + 3.17C/ p , 

e 2 F 2 

where Up — 4m °2 is the so-called ponderomotive poten- 
tial [19, 20J . The difference in the two cases of free graphi- 
nos and active electrons in a graphene layer is mainly due 
to the large excursion-radius in an intense field of the 
active particles which explore very different dispersion 
relations and hence different induced currents. In the 
case considered above, the relative intensity spectrum for 
2D graphene is reminiscent of the high harmonic emis- 
sion spectrum observed recently |13j in experiments with 
bulk ZnO crystals and analyzed within a rather similar 
ID tight-binding model. 



the Fourier transform of the current can be readily ob- 
tained: 



F.T.[J^{ Py ,t)){nw) 



= (-2st&)~(e ifc « 6 - (-l)" e - lfc » h ) 

x(-l)"J„(p) (42) 
= (-2stb) sin (k y b)J n (p) 

for, n=even (43) 
= (— 2stb)icos (k y b) J n (p) 

for, n=odd (44) 



Using the above expression for the current in Eq. ( 34 ), we 



obtain for the relative intensity of the harmonics emitted 
from an initial statefcy in the s th. band analytically in 
the form: 

S x s) (noj) = 4t 2 b 2 (nuj) 2 J 2 (p) sin 2 (kyb), 

for, n=even (45) 
= U 2 b 2 (nuj) 2 J 2 (p) cos 2 (kyb) 

for, n=odd (46) 

It is immediately seen from the formulae above that for 
q y = only odd harmonics are generated while for non- 
zero valucsof q y both odd and even harmonics are emit- 
ted. From the known exponentially decreasing asymp- 
totic behavior of Bessel functions of order n, for n > p, 
one also expects an exponential "cutoff" of the harmonic 
emission above 



""cutoff 



CO 



(47) 



CONCLUSION 



We obtain a class of exact adiabatic solutions of the 
time-dependent Dirac equation for the low energy mass- 
less excitations of 2D graphene sheets or "graphinos" in- 
teracting with an electromagnetic field. They are ex- 
pected to be useful for analyses of non-perturbative pro- 
cesses in intense laser fields. As a first case of their use we 
predict copious emission of harmonic radiation by THz 
laser with considerable relative intensities. Thus, for ex- 
ample, calculation near the Dirac point (e.g. q x a = q y b = 
0.01) shows that many higher order harmonics can be 
generated by a linearly polarized 3 THz laser with a fo- 
cused intensity of 50 MW/cm 2 , whose relative intensities 
may change by only about two orders of magnitude from 
the 3rd to the 81st harmonic order. Furthermore the 
calculated distribution clearly shows two interesting do- 
mains of "revival" and "plateau" formation in the relative 
harmonic intensity. Preliminary results for a laser lin- 
early polarized along the y-axis, and for initial states near 
the Dirac point show that high harmonic generation from 
a 2D graphene layer treated in a tight-binding approxi- 
mation can be qualitatively different from the spectrum 
for the free graphinos. The former, instead, is rather 
analogous to the recently observed high harmonic spec- 
trum from bulk ZO crystals. 
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